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Motivated by the work of Kitaev, we construct an exactly soluble spin-| model on honeycomb lattice whose 
ground states are identical to A\ x p x + AiyPy + i(A2 X p x + A2j,p H )-wave paired fermions on square lattice, 
with tunable paring order parameters. We derive a universal phase diagram for this general p-wave theory 
which contains a gapped A phase and a topologically non-trivial B phase. We show that the gapless condition 
in the B phase is governed by a generalized inversion (G-inversion) symmetry under p x <-> -^-Py- The 
G-inversion symmetric gapless B phase near the phase boundaries is described by 1+1 -dimensional gapless 
Majorana fermions in the asymptotic long wave length limit, i.e. the c = 1/2 conformal field theory. The 
gapped B phase has G-inversion symmetry breaking and is the weak pairing phase described by the Moore- 
Read Pfaffian. We show that in the gapped B phase, vortex pair excitations are separated from the ground state 
by a finite energy gap. 



I. INTRODUCTION 

The low energy excitations of a topologically non-trivial 
phase have remarkable properties. A well-known example 
is the quasihole excitation of the Laughlin state in the frac- 
tional quantum Hall effect (FQHE), which carries fractional 
charge and anyon statistics. The most intriguing possibil- 
ity of a topological phase of matter is the nonabelian FQHE 
proposed by Moore and Read 1 for even-denominator filling 
factors, e.g. v — | 2 . The quasiparticle excitations, vor- 
tices in the Moore-Read Pfaffian wave function, have non- 
abelian statistics 1 which plays a fundamental role in topologi- 
cal quantum computation 3,4 . The key ingredient in the Moore- 
Read Pfaffian state is that the topologically nontrivial part of 
the wave function is asymptomatic ally the same as the pair 
wave function in a p x + ip y -wa\e fermion paired state 6 in the 
weak pairing phase 7 . The existence of the exotic non-abelian 
statistics is thus likely a generic property of the more tangible 
time-reversal symmetry (T-symmetry) breaking p-wave pair- 
ing states. 

Recently, Kitaev constructed a spin-^ model with link- 
dependent Ising couplings on the honeycomb lattice 4 -. Ki- 
taev showed that the model is equivalent to a bilinear Majo- 
rana fermion model and is thus exactly soluble. A topologi- 
cal non-trivial gapless phase (the B phase) was discovered. (A 
Jordan- Wigner transformation to a model with two-Majorana 
fermions for this model has been proposed in 5 ). In the pres- 
ence of a T-symmetry breaking term, the B phase becomes 
gapped and exhibits vortex excitations obeying nonabelian 
statistics. The model also has a topologically trivial, gapped 
A phase. The two phases are separated by a topological phase 
transition via a gapless critical state. These properties strongly 
resemble the weak and strong pairing phases and the critical 
state in the p x + ip y -wave paired states of spinless fermions 7 . 

The interconnections among the Kitaev model, the p-wave 
paired fermions, and the Moore-Read Pfaffian and its excita- 
tions have not been well understood previously. In particular, 
it is important to understand the universal properties among 
these systems, analogous to finding the universality class in 
statistical mechanics models. In this paper, we show that 



the Kitaev model is a special case of a broader class of two- 
dimensional spin-i models whose ground states are equiva- 
lent to general paired fermion states in the p-wave channel. 
Indeed, the vortex-free Kitaev Hamiltonian maps to an ex- 
act BCS fermion pairing model with i(p x + p y )-wave attrac- 
tions on a square lattice 9 . Our generalized model includes 
both the p x + ip y wave paired states and the original Ki- 
taev model as special limits. It is an exactly soluble model 
with minimal three and four-spin interactions. We show that 
the vortex-free ground states of this model are described by 
A lx p x + A ly p y + i(A 2x p x + A 2y p y )-wave paired fermion 
states with tunable pairing order parameters A a j, on a square 
lattice. We find that the structure of the phase diagram is de- 
termined by the geometry of the underlying Fermi surface. 
It contains both topologically trivial (A) and nontrivial (B) 
phases. The A phase is always gapped and corresponds to 
the strong pairing phase. The B phase can be either gapped 
or gapless even if T-symmetry is broken. We find that gap- 
less excitations in the B phase is protected by a generalized 
inversion (G-inversion) symmetry under p x <-> ^f-Py and the 
emergence of a gapped B phase is thus tied to G-inversion 
symmetry breaking. For instance, the p x + ip y wave paired 
state is gapped while p y + zp y -wave paired state is gapless 
although they both break the T-symmetry. The critical states 
of the A-B phase transition remains gapless whether or not 
T- and G-inversion symmetries are broken, indicative of its 
topological nature. Indeed, if all A a b are tuned to zero, the 
topological A-B phase transition is from a band insulator to a 
free Fermi gas. The Fermi surface shrinks to a point zero at 
criticality. 

We show that the gapped B phase is a weak pairing state 
while the G-inversion symmetric ground states are extended. 
The gapless phase was not well-understood before. We show 
that the effective theory near the phase boundary corresponds 
to 1+1 -dimensional massless Majorana fermions in the long 
wave length limit, i.e., a c = 1/2 conformal field theory or 
the 2-dimensional Ising model. The vortex excitations are im- 
portant in the family of Kitaev models since the vortex exci- 
tations may obey anyon statistics 4 . The vortex excitation en- 
ergies have been numerically estimated in the A phased and 
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the B phase 11 . We study the vortex excitations in the gapped 
B phase in the continuum limit and show that the vortex pair 
excitations cost a finite energy. This is consistent with the 
results of numerical calculations 11 and suggests that vortex 
excitations may have well-defined statistics. 



n. THE GENERALIZATION OF KITAEV MODEL 

We extend the Kitaev model on the honeycomb lattice by 
introducing minimal three- and four-spin terms in the Hamil- 
tonian, 
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where <j x,y ' z are Pauli matrices, x-,y-,z-links are shown in 
Fig. [Tapper panel), 'w' and 'b' label the white and black 
sites of lattice, and e x ,e yi e z are the positive unit vectors, 
which are defined as, e.g., ei2 = e z ,e23 = e x ,eei = e y . 



and \ x , y are tunable real parameters. The orig- 



inal Kitaev model has n a = X a = 0, a = x, y. Adding a 
T-symmetry breaking external magnetic field corresponds to 
k x = K y = k ^ and a K 2 -termi£. It is important to note 
that the generalized Hamiltonian maintains the Zi gauge sym- 
metry acted by a group element, e.g., 
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with [H, Wp] = 0. In fact, one can construct Zi gauge 
invariant spin models with higher multi-spin terms, e.g., 

cr 9 f7 10 cr l cr 2 cr 3: cr 9 cr 10 cr l cr 2 o '3 cr 4 and °9 cr 10 cr l °2 and 

so on. One can also add the 'z'-partners of K x<y and \ XtV 
terms so that the model becomes more symmetric. However, 
adding these term or not will not affect our result in this paper 
as we will explain later. 

We now write down the Majorana fermion representation of 
this spin model. Let b x y z and c be the four kinds of Majorana 
fermions with b 2 x z = 1 and c 2 = 1. The spin operator is 
given by 

i 1 

o- a = 2^ qC_ 2 e abcbbb c ). 

Restricting to the physical Hilbert space, one needs to require 4 




FIG. 1: (Color online) Upper panel: The honeycomb lattices and 
links . Lower panel: The vortex excitations. The empty circles de- 
note Wp = 1 and the filled circles denote vortices with Wp = — 1. 



and thus a a = ib a c. The Hamiltonian now reads 
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Where K^ h+Sz — Kx^He^Hez+e^He,' ^b,b+2e z +e m ~ 
^xK,b+e z U b+e z ,b+e z +e lc uZ b+e z +e :L ,b+e z +e :c +e z etC and U ij = 

ibfbl on a-links. It can be shown that the Hamiltonian com- 
mutes with ufj and thus the eigenvalues of ufj — ±1 because 
(ufj) 2 = 1. Since the four spin terms we introduced are 
related to the hopping between the 'b' and 'w' sites, Lieb's 
theorem 8 is still applicable. The third spin terms are 'b' to 
'b' and 'w' to 'w' and Lieb's theorem is not directly ap- 
plicable. However, according to Kitaev^, one can still take 
u bw = — ii° b = 1. The Hamiltonian for the ground state free 
of the Zi vortices (Wp — 1 for all P) is given by 
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where s represents the position of a z-link, A Q — — o 
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and 
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where E p = y£,p + (Ai, p ) 2 + (A2, p ) 2 is the dispersion, 
A p = Ai.p + «A2. P , and (u p ,v p ) are the coherence fac- 
tors with |u p | 2 = |(1 + ^),|w p | 2 = 1(1 - ^) and 

V"p = -(£ P -£p)/Ap. 



FIG. 2: Left panel: The effective square lattice. Right panel: Phase 
diagram in the space spanned by (J x , J y , J z ). This is a (1,1,1) cross 
section with all J's positive. 



III. MAPPING TO A p-WAVE PAIRED STATE 

The next step is to map the Majorana fermion Hq to a 
fermion model. Defining fermions on the z-links by 9 

d s = (c s ,6 + ic StW )/2, d\ = (c s , b - ic SjW )/2, 

H n becomes 

Ho = Jz £(djd. - 1/2) + JMde+e, - dJ s+e J 

s 

+ \ x ^2(dl +e jt - d s+e J s ) (4) 

S 

+ iKx + d\dl +e J + y partners. 

s 

This is a quadratic model of spinless fermions d s on the square 
lattice (Fig. [2] left panel) with general p-wave pairing. If we 
include the V partner of the three and four spin terms in eq. 
(H), we have additional corresponding terms in eq. which 
are the next nearest neighbor terms in the square lattice. These 
terms will not qualitatively affect our result. Returning to Ma- 
jorana fermions, a link fermion is a superposition of two Ma- 
jorana fermions in a link. Therefore, the paring of the link 
fermions reflects the 'pairing' of Majorana fermions. 
After a Fourier transformation, eq. (|4]i becomes 

#o = j2^ d l d p + ^ R ( d l d - P + d P d -p) 
p 

+ i^(dldl p -d p d- p ) (5) 
where the dispersion and the pairing functions are 

£p = J z - Jx cosp x - J y cospy, 

A QiP = A ax sin ^ + A ay sinp y , a =1,2 

with &i, x ( y ) = n x ( y) and A 2)X (y) = K( y )- We have thus 
shown that the ground state of the extended Kitaev model in 
Eq. (1) are equivalent to general p-wave paired fermion states. 
The quasiparticle excitations are governed by the BdG equa- 
tions 

EpU p = £pu p - A*Up, E p v p = -£pv p - ApMp (6) 



IV. PHASE DIAGRAM IN TERMS OF THE p-WAVE 
STATES 

We now turn to the properties of this p-wave paired state 
of link fermions. It is instructive to consider the free fermion 
dispersion. The condition £ p = defines a topological transi- 
tion between a band insulator and a metal with a Fermi surface 
in the absence of pairing. The solution is given by | cosp* | = 
cosp;| = 1, p* = (0, 0), (0, ±tt), (±tt, 0), (±tt, ±tt), where 
J z ± J x ± J y — 0. Without loss of the generality, one con- 
siders only J x y > and J z > 0. Then, £ p » = corresponds 
to the inner triangle of the (l,l,l)-cross section in ( J x , J y , J z ) 
space (see Fig. 2 (right panel)). Notice that the p-wave pairing 
gap functions A a p vanish at p* and therefore, this triangle is 
the gapless critical boundary separating the A and B phases. 
Outside the triangle, £ p > 0. Thus the A phase is gapped. In 
the limit p — » p*, the pair correlation g p = v p /u p is analytic 
near p*, implying tightly bound pairs in positional space and 
hence the A-phase as the strong-pairing phase (See below for 
detailed discussions) 7 . The global structure of the phase dia- 
gram is invariant in the generalized (J x , J y , J z ) space because 
our minimal three- and four-spin extension does not change 
the topology of the underlying Fermi surface. 

The nature of the B phase is much more intriguing. In- 
side the triangle, £ p , Ai p and A2. p can be zero individ- 
ually. The gapless condition (E p = 0) requires all three 
to be zero at a common p*. This can only be achieved if 
(i) one of the A a p = or (ii) Ai, p oc A2 ;P . If either 
(i) or (ii) is true, ^ p and A p can vanish simultaneously, i.e. 
E p = at p*, and the paired state is gapless. Otherwise, the 
B phase is gapped. Note that contrary to conventional wis- 
dom, T-symmetry breaking alone does not guarantee a gap 
opening in the B phase. The symmetry reason behind the gap- 
less condition of the B phase becomes clear in the continuum 
limit where E p = implies that the vortex-free Hamiltonian 
must be invariant, up to a constant, under the transformation 

Px <-* VPy an< l Jx V 2 Jy w i tn V = 3^ w i tn a = 1 or 
2 and for nonzero A. We refer to this as a generalized in- 
version (G-inversion) symmetry since it reduces to the usual 
mirror reflection when r] = 1. This (projective) symmetry 
protects the gapless nature of fermionic excitations and may 
be associated with the underlying quantum order—. Kitaev's 
original model has Ai,j = 0, and is thus G-inversion invari- 
ant and gapless. The magnetic field perturbation 4 breaks this 
G-inversion symmetry and the fermionic excitation becomes 
gapped. A special case with G-inversion symmetry breaking 
is A p oc sin pa; + i smp y , i.e. the p x + ip a -wave paired state 
discussed by Read and Green in the continuum limit 7 . 

The iV-fermion ground state wave func- 
tion in the general p-wave paired state can 
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be written down as a Pfaffian for N even: 
¥(!•!, ...,r x ) =l/(2 w / 2 (7V/2)!)E P sgnPnt / i 2 5(rp 2 ,_ 1 - 
rp 2i ) where g(r) is the "pair correlation", i.e. the Fourier 
transform of g p — v p /u p in the BCS wave function 

l«> = n p l«pr /2 exp(iEp5p4^p)l°)- The wave 
function exhibits very different behaviors in the long wave 
length limit in different parameters 13 . In the A phase, £ p > 
as p — > , thus <7 p oc A p . The analyticity of g p leads to 
<7(r) oc e~ Air as in the strong pairing phase of a pure p x + ip y 
statei In the gapped B phase with G-inversion symmetry 
breaking, £ p < as p — > 0. Defining p[ = A a iPi with 
a = 1, 2 and £ = a;, y, it follows that g p oc —, , leading 

to <?(r) = p with x' a — A~^Xi and thus a weak-pairing 

phase. Identifying z' = x[ + ix' 2 , we see that the ground 
state of the gapped B phase corresponds exactly to the 
Moore-Read Pfaffian. It is easy to show that, u and v* obey 
the same BdG equation in this general p-wave paired state, 
such that the anti-particle of the quasiparticle tp = (u, v) is 
itself, i.e., it is Majorana fermion obeying Dirac equations in 
2+ 1 -dimensions^. 

We now discuss the nature of the gapless B phase in the 
general model with G-inversion symmetry. In this case, E p = 
at p = ±p* which are the solutions of £ p = and, say, 
A p = Ai.p = 0. At p*, the fermion dispersions are gen- 
erally given by 2D Dirac cones. However, by a continuous 
variation of the parameters, one can realize a dimensional 
reduction near the phase boundary where the effective the- 
ory is in fact a 1+1 dimensional conformal field theory in 
the long wave length limit. Let us consider parameters that 
are close to the critical line with | sinp* <c | cosp*| where 
,9 q = sgn[q x A lx cosp* + q y A ly cosp*} = sga(q' x ) with 
q = p p*. Doing the Fourier transform, we find 

9 (t) = J d q ' x dq' y e^' + <y S gn( q ' x ) 

= S(y') fd q ' x ^ S m q > x x'~ 5 -^±. (7) 

J \1x I x 

The (5(y)-function indicates that the pairing in the gapless B 
phase has a one-dimensional character and the ground state is 
a one-dimensional Moore-Read Pfaffian. The BdG equations 
reduce to 

id t u = -£Ai x (l + irj)d X 'V, id t v = £Aia,(l - ir])d x >u, (8) 

with r\ = -^ L . Thus, the gapless Bogoliubov quasiparti- 
cles are one-dimensional Majorana fermions. The long wave 
length effective theory for the gapless B phase near the phase 
boundary is therefore the massless Majorana fermion theory 
in 1+1 -dimensional space-time, i.e. ac — 1/2 conformal field 
theory or equivalently a two-dimensional Ising model. 

V. TOPOLOGICAL INVARIANT IN A AND B PHASES 

We note that there is no spontaneous breaking of a contin- 
uous symmetry associated with the phase transition from A 
to B phases. Kitaev has shown that the A phase in his model 



is topologically trivial and has zero spectral Chern number, 
while the gapped B phase is characterized by the Chern num- 
ber ±1—. This fact was already discussed by Read and Green 
in the context of p x + ip y paired state. Here we follow Read 
and Green 7 to study the topological invariant in a general p- 
wave state. 

In continuum limit, p = (p x ,p y ) lives in an Euclidean 
space R 2 . However, the constraint |u p | 2 + \v p \ 2 = 1 pa- 
rameterizes a sphere S 2 . As |p| — > oo, £ p — * E p such that 
Vp — » 0. Therefore, we can compactify R 2 into an S 2 by 
adding oo to R 2 where v p — ► 0. The sphere |m p | 2 + |t> p | 2 = 1 
can also be parameterized by a pseudospin vector n p = 
(Ai, p , -A 2 ,p,£p)/-E p because |n p | = 1. (u p ,v p ) thus de- 
scribes a mapping from S 2 (p € R 2 ) to S 2 (spinor |n p | = 1). 
The winding number of the mapping is a topological invari- 
ant. The north pole is u p = l,v p = at |p| = co and 
the south pole is u p — 0, v p = 1 at p = 0. In the n p 
parametrization, no = (0, 0, ^jr) = (0, 0, 1) at |p| = co and 

(0, 0, %) = (0, 0, —1) at p = 0, corresponding to either the 
north pole or south pole. 

In the strong pairing phase, we know that u v — > 1 and 
v p — > as p — > (or equivalently, £ p > 0). This means 
that for arbitrary p, (u p , v p ) maps the p-sphere to the upper 
hemisphere and the winding number is zero. That is, the topo- 
logical number v = in the strong pairing phase. 

In the weak paring phase, u p — » and v p — > 1 as p — > 
0. This means that the winding number is nonzero (at least 
wrapping once). For our case, the winding number can be 
directly calculated and is given by 

v = 4^ / dp x dp y n p ■ (<9 Px n p x d Py n p ) = 1 (9) 

Defining P(p) = h(l + n p • a), which is the Fourier compo- 
nent of the projection operator to the negative spectral space 
of the Hamiltonian, this winding number can be identified as 
the spectral Chern number defined by Kitaev 4 

v = ^-\ T4P-(d Px P-d Py P- - a Pu P-a p .P-)]d Ps dpQ0) 

where P_ = I — P is a projective operator. This spectral 
Chern number vanishes in the strong pairing A phase but takes 
an integer value in the weak paring B phase. Thus, the phase 
transition from A to B is a topological phase transition. 

VI. VORTEX EXCITATIONS 

We have discussed the universal behaviors of the ground 
state. We now turn to discuss the Z 2 vortex excitation in 
the spin model which corresponds to setting Wp = — 1 for 
a given plaquette. The Hamiltonian in the Majorana fermion 
representation is bilinear and the energies of the vortices can 
be estimated both in the A phase 10 and the B phase 11 . How- 
ever, it remained difficult to obtain analytical solutions of the 
wave functions with two well-separated vortices. We have 
shown that the ground state sector is equivalent to the p x + ip y 
pairing theory for fermions on the square lattice. Therefore, 
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the Pfaffian state is the ground state wave function in the con- 
tinuum limit in the weak pairing phase. Our strategy is to 
evaluate the energy of the trial wave function containing vor- 
tices above the Pfaffian state in the continuum limit. For two 
well separated half-vortices located at w± and w 2 shown in 
Fig. 1 (lower panel), the Moore-Read trial wave function has 
been well-studied 1 ' 7 and is given by 

...z N ; iOi,toa ) cx Pf(g'(zi, Zj;wi,w 2 )), (11) 
(zi - wi)(z 2 - w 2 ) + (wi <-> w 2 ) 



f(zi,z 2 ;wi,w 2 ) oc 



zi - z 2 



The second quantized state corresponding to this wave func- 
tion reads 

\w!,w 2 ) oc exp{i ^2 9'(zi,Z2;w 1 ,w 2 )dl 1 dl 2 } (12) 



Performing a Fourier transformation, we have 

\wi,w 2 ) ocexp{-^^(K)4c+k4c-k}' 

K,k 

where k = ki — k 2 and K = ki + k 2 are the relative and the 
total momenta of the pairs and g k (K) is the Fourier transform 
of g'(rx,r 2 ). One can show that, 



«i(K) ~ ( T _ 
4i 



A 



■)6<K) 



1 

k wit«2|k| 2 fc ; 
B (w 1 +w 2 )C, 



k y Wl w 2 \K\ 2 K 2 wiw 2 \K\ 2 K' 
= g' k S(K) + lg(K) 



(13) 



where A, B and C are positive constants and g(K) is inde- 
pendent of k. Thus, 

\w!,w 2 ) ocexp{i^^4 d -k+^( 1 / fc )3'( K ) d K+k rf K-k} 

k K,k 

Such a vortex pair is shown in Fig. 1 where the red z-links 
have Ubw = — 1 and all others have Ub w = 1. The correspond- 
ing Hamiltonian can be written as H = Ho + 5H, where Ho 
is the vortex-free Hamiltonian and 5H is the vortex part. The 
latter is expressed as a sum of the pairing and chemical poten- 
tial terms according to Eq. (5) over the red z-links extending 
in the ^-direction (i.e., the direction with x = y) between the 
vortices. It is straightforward to show that 5H has the follow- 
ing expectation value in the vortex state, 



(wx,w 2 \5H\wi,w 2 ) 



(14) 



JW2 (P£+Pj) x 



P(<Pe 
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where /(p{,Pj) is an analytical function of p^ + p^. This 
means that there are no a continuum spectrum above the vor- 
tex pairs and then the vortex pairs are also separated from 
other higher energy excitations. On the other hand, one can 



check that since [H , £ Kk 9k( K f °)4c+k4c-k] = °> the 
K^O sector does not play a nontrivial role in calculating the 
energy E v of such a vortex pair. The latter is given by 

E v = (wi,w 2 \H\wi,w 2 ) = (ioi,w 2 |i?o|wi,W2) 
= E k \u k Sg k \ 2 (wiw 2 \dudl\wiw 2 ) 



Y,E k \u k dg k \y(l+\g k r), 



B 



k W1W2K 2 



where sj? = g' k + \{ 
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— i/k- i " 

KT) = V(i 



9 k — 5k- Physically, the factor {wiw 2 \d k d\\w\W 2 ) 



I5f| 2 /(1 



(15) 

)|k->o and Sg k 
1 

\g' k \ 2 ) in Eq.(15) is the quasi- 



(w 1 +w 2 )C 



hole distribution when the two vortices are located at w\ and 
w 2 . Thus, E v indeed corresponds to the energy cost to ex- 
cite the vortex pair. We have evaluated the vortex pair en- 
ergy E v in different limits. First, if 101 and w 2 were sent 
to infinity before K — > 0, then g'P —> g k and we recover 
the ground state. Second, if K — > while w\ and w 2 re- 
main finite, then E k \u k 5g k \ 2 / {1 + \g' k \ 2 ) ~ E k \u k \ 2 , which 
tends to |k| 2 in both the small and large k limits. The exci- 
tation energy is thus high but finite due to the short distance 
cut-off. (Recall that £ k < in the B phase). It is indepen- 
dent of wi and w 2 and as a result the vortices are decon- 
fined. The third case is when the vortices are far from the 
origin such that |-K"ioi| and |-ftTio 2 | are finite. In the short dis- 
tance, large k limit, E k \u k 5g k \ 2 /{I + \g' k °\ 2 ) ~ |k|~ 4 . On 
the other hand, in the long wavelength limit with |k| — > 0, 
\u k \ ~ |k|, \g' k °\ — > |k|~ 3 and E k — > constant such that 
E k \u k 5g k \ 2 /{1 + \g' k °\ 2 ) ~ |k| 2 -> 0. As a result, the vor- 
tex pair energy E v is free of infrared divergences and is only 
weakly dependent on w\ — w 2 . Therefore, the vortices are 
also deconfined. Finally, since E k increases with the paring 
parameters n and A in Eq.(4), E v is expected to increase with 
increasing paring gap parameters. These results are consistent 
with the finite size numerical calculations of the vortex pair 
energy in the lattice model 11 . Our analytical results suggest 
that the vortex pair described by Eq.(l 1), while costing a high 
energy in the bulk, corresponds to low energy excitations near 
the edge of the system. 

We note an important difference between this p-wave the- 
ory and a conventional p-wave superfluid: Instead of spon- 
taneously breaking the U(l) symmetry in an usual p-wave 
superfluid, only the discrete Z 2 symmetry is broken and the 
U(l) symmetry is absent in the present model. The vortices 
studied here are thus Z 2 vortices instead of U(l) vortices. As 
a consequence, in the gapped B phase, the vortices are in the 
deconfinement phase 15 instead of being logarithmically con- 
fined in the p-wave superfluid. This fact can be easily seen be- 
cause &-ax0y) are rea l an d there is no U(l) phase factor whose 
gradient gives rise to a vector field of the vortex. Our analysis 
of the vortex pair energy also shows this difference. 

The finiteness of E v and the vanishing of (SH) = im- 
ply that the vortex excitations are separated either from the 
ground state or other excitations. This is consistent with anal- 
ysis of Read and Green on the U(l) vortex excitations in the 
p-wave paired state 714 . To determine how close the Moore- 
Read vortex state is to the exact vortex excitations in this 
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model requires a numerical calculation of the overlapping be- 
tween the exact eigenstates and the Moore-Read vortex wave 
functions. A more important question is the realization of the 
Read-Moore four-vortex state which has a two-fold degener- 
acy with the vortices obeying non-abelian statistics 1 . We leave 
these studies to future works. 



VII. CONCLUSIONS 

We have constructed an exactly soluble spin model with 
two-, three- and four-spin couplings on a honeycomb lattice. 
The ground state sector of this model on the honeycomb lat- 
tice is mapped to a p-wave paired state of the link fermions 



on a square lattice with general pairing parameters. Based 
on the general p-wave paired states, we analyzed the phase 
diagram of the system and the properties of topologically dif- 
ferent phases. We found that our phase diagram is universal 
and includes both the Kitaev model and the Pfaffian state in 
its universality class. 
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